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I. INTRODUCTION 



The most important feature of the matter -antimatte r interaction for 
rocket propulsion is matter -antimatter annihilation. The occurrence of 
annihilation is necessary for the release of energy for propulsion, but the 
avoidance of annihilation is necessary for antimatter storage. In both cases 
a knowledge of annihilation rates under various circumstances is important. 

The principal purpose of this report is to consider how the annihilation rates 
are affected by the various circumstances and how this, in turn, affects the 
choice of the means of antimatter storage- 

Only a very limited number of accurate calculations that may be used to 
determine accurate annihilation rates have been carried out, and it is beyond 
the scope of this study to carry out such calculations. Therefore, in consid- 
ering the feasibility of various means of antimatter storage, upper and lower 
limits for annihilation rates were employed. For a given means of storage, a 
sufficiently high lower limit on the annihilation rate renders the means unfea- 
sible, or a sufficiently low upper limit on the annihilation rate renders the 
means feasible. 

To determine limits on the annihilation rates, upper and lower bounds tn 
atom-antiatom interatomic potential energies and rearrangement-annihilation 
cross sections were obtained. The upper bounds to the interatomic potential 
energies were obtained by use of a simple atomic model that was also used to 
obtain bounds on other important quantities. The lower bounds to the rearrange 
ment cross sections turned out to be so low that other means could be employed 
to obtain lower bounds to the annihilation rates. Consequently, the lower 
bounds for the cross sections were not employed explicitly in the determination 
of annihilation rates for various means of storage. 

Section III of this report contains a general description of the processes 
that are involved in matter- antimatter annihilation. In Section IV, the impor- 
tant features of atom- antiatom interatomic potential energies for various types 
of atoms and antiatoms are described, and a model is developed that allows 
the approximate calculation of the upper bound to the interatomic potential 
energy for any atom-antiatom pair. The lower bounds are also obtained. 

Based on these potential energies, formulae for the upper and lower bounds 
for atom -antiatom annihilation cross sections are obtained in Section V. 
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Equations for the annihilation energy production rate and the annihilation 
degradation rate for matter -antimatter solids in contact and mixed solids are 
obtained in Section VI by using the model of Section IV. These equations are 
applied in Section VII in considering the possible existence of kinds of matter 
and antimatter that could exist as a stable, solid, homogeneous mixture, and 
are applied in Section VIII to the problem of solid matter and antimatter in 
surface contact. Lower bounds to annihilation energy production rates for a 
gaseous matter -antimatter mix are obtained in Section IX; and in Section X, 
the annihilation problems associated with high vacuum, no- contact, antimatter 
storage are considered. In Section XI, suggestions for future studies are 
made. 

It was not the purpose of this study to consider means for the production 
of the requisite amounts of antimatter. Consequently, no specific composition 
or physical state was assumed for the antimatter. 

Atomic units, in which e = a Q = h = m e = 1,. are used in this report, 
except where stated otherwise. 

II. LITERATURE AND OTHER INVESTIGATIONS 

Past and present research on the aspects of the matter-antimatter 
interaction that are relevant to rocket propulsion is very limited. 

The most important publications are as follows: 

A paper by Morgan and Hughes (Ref. 1 ) that gives detailed general 
and specific consideration to most important atom-antiatom 
interactions at thermal energies. 

An earlier paper by Morgan and Hughes (Ref. 2) that includes 
consideration of the matter-antimatter interaction at very high 
energies. 

A paper by Junker and Bardsley (Ref. 3) that presents the best 
calculations to date on the thermal energy hydrogen-atom (H) 
antihydrogen-atom (H) interaction. 

A paper by Omidvar (Ref. 4 ) that gives the results of a Born 
approximation calculation for the H-H rearrangement and annihi- 
lation cross section that is important for kinetic energies some- 
what above the normal thermal energy range. 
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The past work of G. Steigman on the positronium formation cross section 
in an H-H collision at energies above normal thermal energies appears in his 
dissertation (Ref. 5). The work of H. Alfven and O. Klein (Ref. 6), Nauenberg 
and Ruderman (Ref, 7), and J. L. Puget (Ref. 8) as it applies to this study is 
superseded by References 1 through 4. 

Positron scattering, which is a subject of considerable theoretical (Ref. 9) 
and experimental investigation, may be of some relevance in future studies to 
the use of antimatter annihilation for rocket propulsion. 

III. ANNIHILATION PROCESSES 

The processes and circumstances that affect the annihilation rate when 
matter and antimatter come into contact occur on three scales: 

(1) Large-scale processes and circumstances that involve the 
composition, physical state, degree of mixing, density, pressure, 
and temperature of the matter and antimatter. Large-scale 
dynamic processes may involve all the complexities of plasma 
physics and thermodynamics, and are not considered in this 
report. Annihilation rates are determined for certain conditions 
of degree of mixing, kinds of atomic species, physical state, 
density, and temperature, usually to see if a given set of condi- 
tions is capable of stationary or nearly- stationary existence. 

(2) Atomic scale processes that have a direct bearing on the probability 
of contact between the individual particles and antiparticles of 
which the matter and antimatter are composed (Ref. 2). Atom- 

. antiatom interactions are of paramount importance here, and their 

consideration forms a large part of this report. Atom-antiatom 
interactions, for general considerations, may be taken to mean 
interactions of neutral or ionized atoms or molecules with neutral 
or ionized antiatoms or antimolecules, although, in specific 
considerations in the following sections, molecules and anti- 
molecules will not be included. 

(3) Nuclear scale interactions wherein the actual, direct particle- 

antiparticle annihilations occur. For cases involving single 

particles or simple nuclei and antinuclei, the annihilation cross 

sections and rates for these processes are fairly accurately known, 

either theoretically or from experiments (Ref. 2). 
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For gaseous matter and antimatter, mixed or in contact, or for 
individually colliding atoms and antiatoms, the types of processes that 
determine the annihilation rates depend principally on the temperature of the 
gases or equivalently on the atom-antiatom collision energy. At temperatures 
and energies in considerable excess of the total atomic binding energy of the 
electrons (e - ) of an atom and the positron (e^*) of an antiatom, the particle- 
antiparticle annihilations occur directly and the annihilation rates are deter- 
mined from the particle- antiparticle annihilation cross sections (Ref. 2). 

This is true irrespective of the state of ionization of the atoms and antiatoms. 
At temperatures and energies around the value of the binding energy down to 
zero the processes affecting annihilation depend on the degree of ionization of 
the atoms and antiatoms. If they are wholly ionized, the particles of the 
nuclei and the antiparticles of the antinuclei annihilate directly. The electrons 
and positrons annihilate directly or indirectly by first undergoing radiative 
capture to form an atom of positronium (bound e+e") and then annihilating. 

When the atoms and antiatoms are not wholly ionized at these temperatures 
and energies, annihilation occurs nearly always in the atoms of nucleonium 
(bound nucleus and antinucleus) and positronium that are formed in atom- 
antiatom rearrangement collisions. Under nearly all circumstances in which 
atoms of nucleonium are formed, the atoms will not subsequently break up 
before annihilation occurs. Because of this and the fact that the rearrangement 
cross sections are considerably greater than the direct particle-antiparticle 
annihilation cross sections, an atom-antiatom annihilation cross section is 
equal to the rearrangement cross section. 

The physics of the atom-antiatom rearrangement reaction is described 
in detail in Ref. 1. Briefly, for each kind of atom-antiatom pair there exists 
a critical radius R c , such that, in a collision, rearrangement will occur if 
and only if the distance of closest approach between the atom and antiatom is 
less than R c . The rearrangement always involves the formation of nucleonium. 
It may not always involve the formation of positronium since the electrons and 
positrons may be emitted as free particles. This latter fact is not important 
for most considerations here where annihilation resulting directly from colli- 
sions is important because the atom-antiatom annihilation energy comes nearly 
entirely from the nucleus-antinucleus annihilation. 
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Consideration of the interaction of solid or liquid matter and antimatter 
is important for solid or liquid surface contact of the two, or for a hypothetical 
homogenous solid or liquid mixture of the two. Here somewhat different anni- 
hilation processes occur. In a solid the atoms and antiatoms are held in fairly 
fixed positions and do not collide. In a liquid the relative positions are not 
fixed, but the distances between atoms and/or antiatoms do not vary greatly 
and are always greater than R c . If a homogenous solid or liquid mixture 
could exist, or if the surfaces of separate pieces of matter and antimatter in 
contact would not interpenetrate, there would have to be a strong atom- antiatom 
repulsion, at least at small interatomic distances. This repulsion would keep 
the nuclei and antinuclei sufficiently far apart to avoid annihilation. However, 
there would always be some degree of overlap between the electron and 
positron wave functions. Annihilation of an atom and antiatom would then 
occur first through direct electron-positron annihilation or by positronium 
formation. Once this occured, the atom and antiatom would have opposite 
charges. The resulting attraction between them or with other neighbors that 
they would polarize would bring them close enough together or close enough 
to others to permit the formation of nucleonium and probably more positronium, 
both of which would subsequently annihilate. 

Because the critical radius R c is defined as the point where two nearly 
parallel exact energy curves join (Ref. 1), it is very difficult to determine. 

It is known to be 0. 639 a Q for an antiproton and a hydrogen atom (Ref. 1) and 
to be slightly less than 1 a Q for a hydrogen atom and an antihydrogen atom 
(Ref. 3). In this study R c was taken to be 1 a Q for all neutral atom-neutral 
antiatom pairs and 0.5 for all atom- antiatom pairs in which one or both of 
the atom and antiatom are ionized. In many and perhaps all cases where the 
atom -antiatom interatomic potential energy is negative for all R, the exact 
value of R c need not be known to calculate atom- antiatom rearrangement- 
annihilation cross sections (Ref. 1). 

IV. ATOM -ANTIATOM INTERATOMIC POTENTIAL ENERGIES 

Consideration is given here and in the following sections to atoms and 
antiatoms that are spherically symmetric and in the ground state. This limits 
the consideration to atoms and antiatoms with one, two, or eight electrons or 
positrons in their outermost shell. The complexities that arise from the 
multipole moments of the charge distributions of nonspherical atoms are too 
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great to consider here. Nevertheless, most of the formulae derived in this 
and the following sections may be applied to such cases, although with less 
accuracy. Magnetic and other relativistic effects, which have only a small 
effect on the interatomic potential energies, are also not considered. 

The atom- antiatom interatomic potential energy, V, is the total energy 
of atom and antiatom with their nuclei held fixed at a distance R apart minus 
the same quantity for R =°o* 

There are two features of V that are present in all atom-antiatom pairs. 
First, there is no interatomic exchange energy as is present in the V of an 
atom- atom pair. , Second, for sufficiently small R, where there is no screening, 
V = -ZZ/R, the coulomb potential energy between the positive nucleus with Z 
protons and, the negative antinucleus with Z antiprotons. 

To describe other features of V it is necessary to distinguish three types 
of atom- antiatom pairs: 

(1) Neutral atom-neutral antiatom 

(2) Neutral atom -antiion (or ion-neutral antiatom) 

(3) Ion-antiion, where the ions are positive and the antiions are 
negative. 

It is helpful to consider a perturbation expansion for V in which the 
perturbation potential energy is the sum of the coulomb energies between the 
particles of the atom and the particles of the antiatom (Ref. 1). In such an 
expansion the lowest (first) order term, Ej, is just the electrostatic potential 
energy between the nucleus and the undistorted electron charge distribution of 
the atom and the antinucleus and the undistorted positron charge distribution 
of the antiatom. The second and higher order terms represent the effects of 
distortions of the charge distributions caused by the other member of the pair. 

For type 3 pairs, E^ and E^ (the second order energy term), and 
probably the sum of the higher order terms, are negative for all values of R. 
Here, V may be described by the equation: 

V = -C(R)/R (1) 

For R 2 2a Q , C(R) is equal to the product of the absolute values of the net 
charges of the ion and antiion. For small values of R, C(R) increases to a 
limiting value of the product of ZZ. There are no accurate calculations avail- 
able for V for type 3 pairs, although the foregoing knowledge of C(R) is 
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sufficient to allow good estimates of ion-antiion rearrangement-annihilation 
cross sections to be made (Ref. 1 ) (see Section V). 

For type 2 pairs, Ej and E^ are again negative for all values of R as is 
V. For R > 4a Q , V is dominated by the monopole -induced dipole energy from 
which is given by 

V.= -a/2R 4 {2) 

where Ot is the polarizability of the neutral member of the pair. Values of a 
for antiatoms are the same as the values of Ot for the corresponding atoms 
which, in most cases, are known. The only accurate calculation of V for a 
type 2 pair is that of Bates (Ref. 10) for p(antiproton)-H (or p-H). For this 
case the long range form of V given by Eq. (2) is sufficient to allow accurate 
approximations of the rearrangement-annihilation cross section to be obtained 
(Ref. 1) (see Section V). 

As discussed in Ref. 1, the character of V for type 2 and 3 pairs, 
principally the fact that it is negative and attractive for all values of R, allows 
simple reliable formulae to be obtained for the rearrangement cross sections 
of these pairs. The formulae provide accurate approximations or at least 
reasonable estimates for the cross sections. Such is not the case for type 1 
pairs. Here Ej is positive and repulsive for R > 2a Q , while (which 
includes the R - ^, R - ®, R - "^, etc. , long range dispersion energy terms ) is 
negative and attractive for all values of R. The relative magnitudes of E^,- 
E 2 * and the higher order terms and their dependence on R determine the sign 
of V as a function of R, but accurate calculations of these quantities are 
difficult for all type 1 pairs. The only accurate calculation of V for a type 1 
pair is the calculation of V for H-H by Junker and Bardsley (Ref. 3). Their 
results indicate that for this pair, at least, the effects of distortion as 
represented by E;? and higher order terms are sufficient to overcome Ej where 
it is positive and render V negative for all values of R, although only barely so 
for R ~ 3a Q . 

Type 1 pairs are of principal importance to considerations involving 
antimatter storage, so it is necessary to obtain at least some knowledge of V 
for pairs of this type for use in the following sections. This may be accom- 
plished through the use of upper and lower bounds to V. For the upper bound, 
Ej may be employed since it is a variational upper bound to V, For the lower 
bound, the long-range van der Waals dispersion energy, which is the leading 
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term in the expansion of for large R, 

V = -c 6 /r6 

along with any reasonable, smooth interpolation between the value of this at 
R ~ 5a 0 and the nucleus -antinucleus potential energy at small values of R may 
be employed. A somewhat better lower bound could be obtained by adding to 
Eq. (3) the R - ® and R - ^ terms, but their presence is of limited value due to 
their relatively small values, and it would greatly complicate subsequent 
calculations. 

Values of for an atom-antiatom pair are equal to those for the 
corresponding atom-atom pair, which are known or may be obtained from 
simple formulae (Ref. 11). The upper bound, Ej, cannot be evaluated as 
easily. is the negative of Ei for the corresponding atom-atom pair, and 

hence in principle could be obtained approximately from a knowledge of 
Hartree-Fock wave functions of neutral atoms. Since the amount of work 
required to do this is too great for a study such as this, a simple atomic 
model was found that allows reasonably accurate values of to be obtained. 

In this model only the outer shell electrons or positrons are considered 
and their wave function is taken to be a product of exponential functions. The 
inner electrons and positrons are placed at the positions of the nucleus and 
antinucleus to form the cores of the atom and antiatom. Hence, the electron 
and positron charge distributions are, respectively, 

Q(r) = - — e and Q(r) = ^ e H (4) 

where r is the distance between the nucleus and a position in the electron 
charge distribution of the atom, and r is similarly defined for the antiatom; 
n and n are the numbers of outer shell electrons and positrons; and p and (3 
are twice the exponent factors in the exponential wave functions. The exponent 
factors in the wave functions are taken to be and J Z Ej , as they are in the 

hydrogenic case, where Ej and Ej are the ionization energies of the atom and 
antiatom. Hence, 

p =^81^ and P = ^8^1 (5) 

This model is equivalent or similar to atomic models that have been 
employed by many authors (particularly in the early years of the quantum 
theory of atoms) to determine the effects on atoms of external fields and other 
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atoms. The accuracy of results obtained from the use of the model depends on 
the specific application and is then usually somewhat uncertain. In this study 
the model was used to obtain reasonably valid bounds to various quantities. 

In this model Ej is the sum of the atom core- antiatom core, atom core-Q, 
Q-antiatom core, and Q-Q potential energies. The result for (3 ^ (3 is 


where 


= nh(0 4- P)/2 e -R' 
1 32R 


(l + ■£) (v - 6v 2 + v 3 ) R 

e vR ^ - (— - 1 ) (v + 4v 2 + v 3 ) 


A 2 3 . 

4v + v ) 


(-^ - l) (v + 6v 2 + v 3 ) 


r’ = & * £ R and v = 


(H ).'(*-» 4) 


and for P = (3 is 


E 1 = 


nnp -R 


R* 


’ [- 1 -tI r ' + tI r ' 2 + ^ R ' 3 ] 


E^ may be expressed as 


_ nn(p + P) » 
E 1 ~ ^2 • E 1 


( 6 ) 


(7) 


( 8 ) 


(9) 


Ej depends only on the variables v and r’ and the characteristics of may 
be readily determined from the characteristics of . 
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1 

The principal characteristics of E ^ are described in terms of the 
quantities depicted in the figure below: 



i t » » i 

The quantities E^ c , E , a , and E are related to the replacement of 

by an equivalent square potential energy barrier. E 1 is nn((3 4- ]3 )/2 times 

E, the energy of a particle incident on this barrier; E* is nn{(3 + (3)/2 

\ _ _ , 

times E , the equivalent height of the barrier; E, is nn((3 + (3}/2 times 

1 c 

E , the height of the inner portion of this barrier; and a is ((3 + ]3)/2 
times a, the width of the barrier. R c = (3 + p/2 R c (R c is the critical 
radius), and the other quantities are self-explanatory. 


The following tables and formulae give these quantities as functions of 
v and the relations between them: 


V 

R 

0 

R 

m 

E 

m 

© 

• 

o 

2.7 

3.74 

9 . 8052 X 10‘ 3 

0 . 2 

2.8 

3.80 

8. 8938 X 10" 3 

0.4 

2 . 9 

3.98 

6.4443 x 10 ‘ 3 

0 . 6 

3.2 

4.36 

3.3037 X 10~ 3 

0.8 

3.9 

5 . 23 

7.6309 X 10‘ 4 
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Nv V 

0 

0. 2 

0.4 

0 . 6 

0.8 

R \ 





0.25 

- 3.3211 

- 3.3072 

- 3 . 2668 

- 3. 2033 

- 3. 1225 

0.4 

- 1 . 8328 

- 1 . 8202 

- 1.7835 

- 1.7260 

- 1.6531 

0.6 

- 1 . 0203 

- 1 . 0098 

- 0. 97905 

- 0 . 93080 

- 0.86826 

1.0 . 

- 0.40620 ' 

- 0 . 39976 

- 0.38084 

- 0 . 35075 

- 0.31237 

1.5 

- 0. 14527 

- 0. 14275 

- 0 . 13504 

- 0. 12186 

- 0. 10363 

2 . 5 

- 0. 0093201 

- 0, 010150 

- 0 . 012114 

- 0. 013724 

- 0. 013027 

4.0 

; 

0 . 0095394 

0. 0087440 

0.0064431 

0 . 0031275 

0.00021887 


\ 

a : 


V 

0 

1 

0.2 

0.4 

0. 6 

i 

0.8 

9. 8 x 10" 3 
8 . 9 x 10" 3 
6. 4 x 10 3 
3. 3 x 10~ 3 

0 

0 

0 

0 


-2 5 

• 1 . 0 x 10 ‘ 

3 . 5 

3 . 8 

3 . 1 

1 . 7 


-3 0 

1.0 x 10 
7.6 x 10" 4 

5.0 

5 . 1 

5 . 2 

5 . 0 

0 

1.0 x 10“ 3 * 5 

6.4 

6.7 

7 . 1 

7.9 

6.3 

1.0 x 10“ 4 *° 

7.9 

8. 2 

9.2 

10.8 . 

12. 2 

-4 5 

1.0 x 10 

9.2 

9. 7 

11. 1 

13.7 

18. 0 

-5 0 

. 1.0 x 10 * 

10. 5 

11. 1 

13.0 

16.6 

23. 7 

-5 5 

1.0 x 10 

11.8 

12.6 

15. 0 

19. 5' 

29.5 

1.0 x 10' 6,0 

-6 5 

13. 1 

14. 1 

17. 0 

22.3 

35.3 

1.0 x 10 

14.4 

15. 5 

19.0 

25. 2 

42. 0 

-7 0 

1.0 x 10 

15.7 

17. 0 

20.9 

27.9 

46.8 
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The quantity E is determined by requiring that the integral of the 
square root of the potential energy of the real and equivalent barriers be 
the same. It is given to a good approximation by the formula, 


eq 


» * 

R - R 
m 


o rp 

' v E m 




B(v) 1 


-l/2(l-v)R^ -l/2(l-v)(a'+R^) 


8 1 - v 


- e 


for v Z 0. 2 


and 


r » i 

R - R , — r 

m o ! — * 


eq 


,/E 


m 


/ , -R^/2 , , -(a'+Ro)/2' 

+ I (r -f 2) e - (a + R + 2} e 

f \ ' m o 

a ' 

for v 5 0. Z 




( 10 ) 


where 


B(v) = - l) (v + 6v 2 + v 3 ) 


For use in Section VII, for purposes of consistency, it is necessary 
to obtain in the form given by Eq. (3) but with calculated by use of 
the above model. This may be done by using the wave functions of this 
model and by using the closure approximation as is done in Refs. 1 and 2. 
The result is 


768 nn 

(n{3 2 + np 2 ) P 2 P 2 


R 


( 11 ) 
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The simple atomic model is also employed in Section Vito calculate 
the e - e + annihilation rate in solid or liquid atomic scale mixtures of 
matter and antimatter. 

V. ATOM-ANTIATOM ANNIHILATION CROSS SECTIONS 


The rearrangement- annihilation cross sections in collisions involving 
the members of type 2 and 3 pairs and an upper bound to the cross sections 
involving the members of a type 1 pair may be obtained by use of the 
semiclas sical method employed in Refs. 1, 2, and 3, which is valid for 
collision energies in the thermal range from about 0. 0G1 eV to a few eV. 

As discussed in Ref. 1, the nature of the lower bound to V for the 
type 1 pair (corresponding to the upper bound to the rearrangement cross 
section) and of V for the type 2 pair is such that only the form of V at large 
values of R is important in determining the rearrangement cross section. 
Hence the forms of V given by Eqs. (2) and (3), which are valid for large 
values of R, may be employed. 


The upper bound to the cross section in the collision of a type 1 pair 


is 


- ■ >• (&) 


1/3 


( 12 ) 


and the cross section in the collision of a type 2 pair is (Ref. 1) 

, 1/2 


-*(¥) 


(13) 


Note that these two results are independent of the value of R c » These 
formulae and Eq. (14) are valid for values of E in the thermal energy range 
(roughly 1 K to 10, 000 K)(Ref. 1). 

For a type 3 pair, V for all values of R is important. Taking V to be 
given by Eq. (1) with C(R) = q q , the product of the charge numbers of the 
ion and antiion, the semiclas sical method yields (Ref. 1) 

<r = (1 + q q/R c 'E) ir (14) 
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for the rearrangement cross section in the collision of a type 3 pair. Due 
to the uncertainty of V for this case and the uncertainty of R c (^0. 5 a Q here), 
Eq. (14) must be regarded as an estimate of the cross section. In Eqs. (12), 
(13), and (14), E is the collision kinetic energy in the center of mass coordin- 
ate system of the pair. 

A lower bound to the rearrangement cross section in the collision of a 
type 1 pair may be obtained by using the upper bound to V, which is E 1 as 
determined from the atomic model calculation in Section IV, and treating 
the scattering as a quantum mechanical barrier penetration problem. 

To do so it is first necessary to determine the effect of the centrifugal 

potential energy barrier on the problem. The centrifugal potential energy is 

given by U = i {£ + 1)/2MR where £ is the angular momentum quantum 

number of the pair and M is the reduced mass of the pair. For small values 

of £ the value of U at R = R c will be less than E^, the maximum value of 

E and will not substantially reduce the amount of the incoming wave that 

reaches R . For large values of £, the value of U at R = R will be greater 
c o 

than and the centrifugal potential energy barrier will substantially 
reduce the penetration. The value of the mean square of the total angular 
momentum L that results in U being equal to at R = R^ is given by 

L 2 = l (£ + 1) = 2ME (15) 

c c m 

It will be assumed the all incoming waves with £ > ^ c no * penetrate to 

R = R c whereas incoming waves with £ < £^ may penetrate to R = R c , 
subject to the effects of the barrier only. If a wave penetrates to R = R c 
it will undergo rearrangement. 

The value of L given by Eq. (15) corresponds to an impact parameter, 

Rj, equal to L/V2ME = \ r E^/ E~ . Hence, in the classical picture, all 

collisions with R^ > \f E^/E ( =R^) do not involve penetration to R c whereas 

collisions with R . < \fE /E will involve penetration to R^ if the B, barrier 
1 rrr c i 

is passed through. The cross section for rearrangement and subsequent 
annihilation is therefore given by 
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TT 


T = 7T 


m 

i 


T 


(16) 


o- = nRp = 


m 


£ 

E 


where T is the coefficient of transmission of the incoming waves through 
the Ej barrier. 

Second, it is necessary to replace the E^ barrier by an equivalent 
square barrier to calculate T. This barrier is shown as the dashed line in 
the figure in the preceding section, and the quantities in terms of which it 
is defined are given there. It is equivalent in the sense that it is the best 
choice of a square barrier to represent the E^ barrier that can be made 
from simple principles. It will not reproduce the effects °f exactly. 

The transmission coefficient T for the square barrier may be deter- 
mined by following the derivation given by Margenau and Murphy (Ref. 12) 
without making the approximation involved in their final result. Doing so, 
the result is, 


T 


4VK 3 



for E < E and 
eq 


T 


4K x /K 3 



(17) 


(18) 
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for E > E where 
eq 


K 

k 

K 

K 


1 


2 


2 


3 


\1 2M(E - E^) 

a = 2a / (p + (3) 

^2M(E eq - E) 

E lc = 2E i c / nn (P + P) 

i 

\l ZM(E~"e ) 

eq 

E = 2E /nn(P + f£) 

eq eq' r ' 

\T 2ME 

E = 2E’/nn(P + F) 


( 19 ) 

✓ 


Equation (16) along with Eqs. (17), (18), (19), and (10) and the tables 

of Section IV constitute a means of obtaining an upper bound for the 

rearrangement- annihilation cross section for a given atom-antiatom pair. 

Employing this means, the lower limit to the H-H rearrangement cross 

• _13 z 

section is found to be 4. 18 x 10 (in units of a = 1 in atomic units) for a 

- 3 ® 2 / 

collision energy of 2 x 10 (room temperature, in units of e /a = 1 in 

atomic units). For the same energy and C^(H-H) = 6.5, Eq. (12) gives 

87.9 for the upper bound, whereas Junker and Bardsley (Ref. 3) find 27.2, 

which is the best value to date. The very low value of the lower bound 

demonstrates the extreme sensitivity of the cross section to the presence 

of a positive maximum in the interatomic potential energy for collision 

energies in the thermal range. The values of the lower bounds to <r for H-H 

and other pairs are so low that in places in the following sections where these 

lower bounds might be employed to give lower bounds on annihilation rates, 

other means were found that gave higher lower bounds on the annihilation 

rates. 


VI. ANNIHILATION RATES IN THE SOLID OR LIQUID STATE 

As was stated in Section III, the first step in atom- antiatom annihila- 
tion in the solid or liquid state is the annihilation of the electrons and 
positrons by direct contact or through positronium formation. Positronium 
formation is difficult to calculate here, and it is energetically impossible 
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for all but a few atoms and antiatoms in a nonscattering situation* In those 
cases it will probably not increase the rate substantially; hence, positronium 
formation is neglected. 

The direct contact annihilation rate may be determined by employing 
the atomic model of Section IV to obtain the total probability for any 
electron of the atom'being at the same position in space as any positron 
of the antiatom when the atom and antiatom are a distant R apart. This 
probability is then multiplied by the appropriate factors (Ref. 13) to obtain 
the total e + - e~ annihilation rate in the atom and antiatom. 

If the result is multiplied by the number of nearest neighbors of the 
opposite character (taken to be 6) and divided by 2 to make it the rate per 
atom or antiatom the result is 



for (3 £ ( 5 ", and 

1. 51 X 10 9 nnp 3 (l + PR + p 2 R 2 /3) e" PR (21) 

-f 

for p = p, where is the annihilation rate in number of e - e annihilations 
alone per second per atom or antiatom. 

The corresponding energy production rates for e + - e" annihilation 
alone are given by 


W = 4.44 x 10 11 ffi 


( 22 ) 
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where W is in calories per mole per second. If the energy from the 
subsequent nucleus -antinucleus annihilation (as discussed in Section III) 
is to be included, W must be multiplied by 183 6 times the mean atomic num- 
ber of the nucleus and antinucleus times the fraction of the nucleus - 
antinucleus mass energy that is released in the nucleus -antinucleus 
annihilation. 

For use in the following section it is important to have a knowledge of 
what R would be in a solid or liquid mixture of atoms and antiatoms. An 

upper bound to R may be found in the following way. The atomic model may 

rsr - 6 
be used to calculate E^ and the long range, R , form for E^. If this is 

done and it is assumed that V = E + E it is indeed found that there is a 

^ ^ - 6 

minimum in V for R equal to several Bohr radii where the R form is a 
good approximation to E^. Thus, if V = E^ + E^ were a good approximation 
to V, solid and liquid mixtures of atoms and antiatoms (neglecting the 
effects of annihilation) could indeed exist at least momentarily with their 
mean atom- antiatom separation being given by this value of R at the "van 
der Walls minimum" of their interatomic potential energy. 

There are three important facts relevant to the accuracy of a V of 
this form. First, in the atomic model it is assumed that p = 2 \T 2E^ and 
p = 2 7 ZEj , where E^ and E^ are the atom and antiatom ionization ener- 
gies. This assumption is expected to be valid when there is only one 
electron or positron in the outer shell, but when there are more than one it 
will underestimate the value of (3 or p by a few percent, as may be shown to 
be true of helium and argon at least. This in turn results in being over- 
estimated. Second, from known values of (Ref, 1 1 ) it may be shown that 
the values of predicted by this model may be incorrect by up to about 
50%. Third, the neglect of perturbation terms in V of higher order than E-, 
almost certainly will yield higher values of V than the exact ones, as the 
calculations of Junker and Bardsley (Ref. 3) indicate for H - H. Taken 
together these facts mean that there may be no van der Walls minimum as 
predicted by the model or, if there is, the model will probably over 
estimate the value of R at which it occurs. 
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Using the model for p = (3 it is found that the value of R for which the 
minimum occurs, gi ven to a good approximation by 


R . = 15. 2 0 

min r 


-0.65 


(23) 


for closed shell atoms and antiatoms and is somewhat less for atoms and 
antiatoms with open shells. For P / P the equation corresponding to Eq. (23) 
has not been derived. However, it may be shown that replacing p by the 
lesser of P and p will result in an overestimate of R . . Hence the equation. 


R . = 15.2 (lesser of p, p) ^ (24) 

will, relative to these and the preceding considerations give an upper bound 

to R . if an R . does indeed exist. (R . does not exist only when V 
min mm ' mm y 

is negative and attractive for all R, In this case R = 0 gives V = - co, so, 

in a sense, R = 0. ) 

mm 

VII. THE STABILITY OF A SOLID OR LIQUID MIXTURE OF 
MATTER AND ANTIMATTER 

An ideal means of storage of antimatter would be in a homogenous 
solid or liquid mixture with antimatter, if a particular mixture could be 
found that was stable over a reasonable time period. If the conditions in 
such a mixture were only somewhat on the side of producing stability, 
annihilation could be commenced by an achievable amount of compression 
or heating. 

Equations (20) through (22) may be used to determine the values of 
R at which the rate of annihilation energy production by e"*~ - e annihilation 
alone is one calorie per mole per second (a value that might conceivably 
be tolerable in storage). The values of R along with lifetimes of the mixtures 
for H - H and helium (He )- antihelium (He) are 
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Mixture 

R 

Lifetime 

H - H 

28. 5 

14, 203 years 

He - He 

22. 0 

14, 203 years 


The above values of R are, however, considerably greater than the 

values of R . . Using Eq. (24) to give R . for substitution into Eqs. (20) 
min min 

through (22), the following results are obtained. For upper limits on 
lower limits on W, and upper limits on the lifetime. 


Lesser 

of 

(3, P 

R . 
min 

W /nn, cal/ 
(mole • s ) 

nn* lifetime 

1 

1 

15, 2 

1. 56 x 10 16 

2.8x10 ^ s 


2 

9.7 

2. 932 x 10 15 

1.5x10 4 s 

(25) 

3 

7.4 

14 

7.734 x 10 

5. 8 x 10" 4 s 



Since all elements have |3' s between 1 and 3, it is therefore clearly 
impossible on this basis to store antimatter in a solid or liquid mixture with 
matter. The addition of the annihilation energy of the nuclei and antinuclei 
would greatly increase the energy production rates. 

VIII. MATTER -ANTIMAT TER SURFACE CONTACT 

Another possible means of storage of antimatter would be to store it 
in contact with solid matter, such as liquid antimatter in a matter bottle, 
or solid antimatter in contact at certain points on its surface with matter. 

To consider this possibility, Eqs. (20) and (21) may be applied to the 
atoms and antiatoms in contact at the surface and may be regarded as 

-f 

degradation rates. One annihilation of an e - e pair is rapidly followed 
by annihilation of the nticleus and antinucleus of the atom and antiatom that 
will take much less time to be completed than the e + - e pair annihilation. 
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Hence, Eqs. (20) and (21) with n and n set equal to 1 give directly the 
number of monatomic atom or antiatom layers that are degraded while 
contact is maintained at the same value of R. It is therefore possible to 
determine the velocity with which the two surfaces must approach the point 
of contact to maintain contact in the face of degradation as well as the rate 
of annihilation energy production per square centimeter. 

Using Eq. (23) for the relation between (3 and the results are 


p 

R . 
mm 

Atomic 
Layers 
Degraded/ s 

Velocity, 

cm/h 

Annihilation Energy, 
cal/ (cm^ * s) 

1 

15. 2 

35165 

10.4 

4.47 xlO 14 (26) 

2 

9*7 

6599 

1. 23 

5. 29 x 10 13 

3 

7.4 

1741 

0. 244 

1. 05 x 10 13 


In calculating the annihilation energy it has been assumed that the 
densities of the matter and the antimatter are one gram per cubic centimeter 
and that the nuclei and antinuclei annihilate completely. In addition the 
formation of a "leidenfrost" layer between the matter and antimatter has 
been neglected. The leidenfrost layer would consist of posit ronium gas 
and gasses of the matter and antimatter, and its presence would lower the 
annihilation rate. If these effects were taken into account and the smallest 
possible area of contact that could support the antimatter were used, it is, 
nevertheless, very doubtful that the rate of energy production could be 
reduced to a tolerable level. It therefore appears impossible to store 
solid or liquid antimatter in contact with solid matter. 

IX. MATTER -ANTIMATTER CONTACT INVOLVING THE 

GASEOUS STATE 

Two possible means of antimatter storage involve the gaseous state. 
The first is a gaseous mixture of matter and antimatter in a magnetic bottle 
or in a matter container. The second is gaseous antimatter in a matter 
container. 
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In the first case we will assume a mixture of equal numbers of atoms 
and antiatoms. If a significant amount of antimatter is to be stored by this 
means, the annihilation occurring on the walls of any container that may be 
present will be insignificant compared to the annihilation within the gaseous 
mixture. At thermal temperatures, the annihilation will come from two 
processes, direct electron positron annihilation as in the solid or liquid 
state followed by rearrangement annihilation involving the resulting ions 
and antiions, and neutral atom-neutral antiatom rearrangement annihilation. 

The atom- antiatom interatomic potential energy is taken to be given 
by V = E^ + E^, where E^ and are determined by the atomic model 
calculation. This gives an upper bound to V for all but large R and a lower 
limit to the annihilation rate. Another consequence of this choice of V is that 
the annihilation due' to the first process will greatly exceed that due to the 
second process. An upper bound to the lifetime of a gaseous mixture can 
then be obtained by making the simplifying assumption that no annihilation 
will take place when the distance between an atom and antiatom less than 
R . - 2a or is greater than R . + 2a and by using Eq. (21) to obtain a 

mean value for .£%. The results may be presented in terms of the ratio of the 
density to gas density at STP required to give a lifetime of one hour to the 
mixture, the diameter of a sphere that would hold one gram of the mixture 
at STP, and the rate of annihilation energy produced per cubic meter 
assuming two nuclei or antinuclei per electron or positron, complete 
nuclei - antinuclei annihilation, and complete trapping of the energy within 
the container: 


p 

p/ p STP 

Diameter, m 

3 

cal/ (meter • s) 

1 

3.4 x.10 " 7 

50 

8 . 97 x 10 4 

2 

4.4 x 10" 6 

11 

9.31 x 10 6 

3 

2 . 9 x 10" 5 

6 

6.26 x 1 0 7 


(27) 


Even if the above annihilation energy rates were reduced by the 
fraction of nucleons annihilated and the fraction of the energy trapped within 
the container or in its walls, they would probably still be prohibitive. If 
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the volumes were increased (along with a proportionate increase in storage 
time) to make the annihilation energy production rates tolerable, the system 
would be so large that it would probably be only barely competitive with 
chemical propellants in terms of the total mass to energy stored ratio. It 
therefore appears that this form of antimatter storage is not practical. 

If the scattering annihilation rates were employed, either with the 
upper limit to the rearrangement cross section for thermal energies 
(Eq. (12)) or with the direct particle- antiparticle cross sections that would 
apply at temperatures of a few thousand Kelvins or greater (Ref. 2), 
then the results would be considerably worse. 

In the second means of storage, annihilation occurs only between the 
atoms of the antimatter gas and the atoms of the container, and the treat- 
ment is otherwise es sentially the same as above. If we assume V = E^ + E^ 
with E^ and E^ being given by the model, then annihilation will occur only 

within a shell about 2a thick which is located a distance of R . from the 

o min 

inside of the container. For various diameters of spherical containers 
and the three values of p, upper limits to the lifetime of the antimatter gas 
and lower limits to the annihilation energy production rate for e - e 
annihilation alone are: ■ 


nn x lifetime (days), energy produced (cal/mole* s)/nn 


n \d 

p\ 

1 cm 

1 m 

10 m ■ 

1 

0. 020, 2. 5 x 10 8 

2. 0, 2. 5 x 10 6 

20. 0, 2. 5 x 10 5 

2 

0. 11 , 4. 7 x 10 7 

11.0, 4. 7 x 10 5 

110 , 4. 7 x 10 4 

3 

0.42 , 1.2 x 10 7 

42,0, 1.2 x 10 5 

420 , 1. 2 x 10 4 


(28) 


These energy production rates are prohibitive and would be more so 
if scattering annihilation had been considered with the cross-section given 
by Eq. (12), or by the high-energy direct particle-antiparticle annihilation 
cross sections (Ref, 2). 
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X. ANNIHILATION ENERGY PRODUCTION RATES FOR 
ANTIMATTER IN A VACUUM 

Another means of possible antimatter storage'is electromagnetic 
suspension of solid antimatter in a vacuum. Annihilation would occur here 
due to vaporization of the antimatter and the surrounding matter structure. 
The composition and positioning of the matter could undoubtedly be chosen 
to make its vaporization insignificant, but the composition of the antimatter 
would, most likely, be determined by other factors. Hence, it is likely 
in a practical situation that significant annihilation would occur only in the 
interaction of vaporized antiatoms with atoms at the surfaces of the 
matter. 

An upper limit to the annihilation and energy production rates may be 
obtained by assuming that each antiatom interacts with the matter surface 
atoms in the same way that it interacts with an individual atom, that 
Eq. (12) for the cross section applies, and that complete nucleus -antinucleus 
annihilation occurs following the rearrangement reactions. This cross 
section and the corresponding interaction are such that for all temperatures 
in the ordinary thermal range, except a few Kelvins or less, several 
angular momentum waves take part in the rearrangement, and the scatter- 
ing is dominated by the inelastic rearrangement process. Under these 
circumstances, the atom-antiatom elastic scattering cross section is about 
equal to the rearrangement cross section (Ref. 1). Thus, when an antiatom 
strikes a matter surface, it will undergo annihilation with a surface atom 
with a probability of about one half. Since the elastic scattering will be 
predominantly in a forward direction, most of the elastically scattered 
antiatoms will enter the matter and subsequently annihilate. It may 
therefore be assumed that all antiatoms striking a matter surface will be 
annihilated. 

It follows that there is a very simple relation between the lifetime 
against evaporation in years, t, of the antimatter, and the upper bound to the 
annihilation energy production rate, W. If it is assumed that the fraction 
of the annihilation energy that is deposited in the antimatter is one half (an 
upper limit since roughly one half of the energy goes into neutrinos and 
most of the energy could be deposited elsewhere) the relation is: 
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(29) 


W = 6. 83 x 10 5 


cal 

mole • s 


where A is the atomic number of the antimatter. Considering the nature of 
W as an upper bound and small values of A, this form of storage would be 
practical [W5 1 cal/ (mole* s) ] if the antimatter has a lifetime against 
evaporation of about one million years or greater, and may be practical 
for smaller lifetimes. Such lifetimes can be readily achieved with matter. 
They could, in principle, be achieved with antimatter; but whether they 
could be achieved or not, in a practical situation, is difficult to determine. 

XI. SUGGESTIONS FOR FURTHER STUDY 


The results of this study were based on the interactions of spherically 
symmetric atoms and antiatoms that possess no electric or significant 
magnetic multiple moments. A better knowledge of the interatomic 
potential energy for such species would lower the gap between the upper 
and lower bounds determined for various quantities. This would be par- 
ticularly important in evaluating the possible means of transferral of the 
antimatter to the combustion chamber of a rocket and its annihilation in the 
combustion chamber. 

The interactions relevant to annihilation of atoms and molecules that 
possess electric or magnetic multipole moments may differ significantly 
from the interactions considered in this study. Multipole moments may 
result in different annihilation rates, and they may conceivably alter the 
atom-antiatom interactomic potential energy to the extent that the mean 
atom- antiatom separation in a solid, atomic scale matter -antimatter 
mixture is sufficiently great to make such a mixture stable over a long 
period of time. * 

Electric and magnetic multipole moments are possessed by some 
atoms and molecules. For these or other atoms and molecules moments 
could be enhanced or induced by the application of external electric and/or 
magnetic fields. The effects of multipole moments on annihilation rates and 
atom- antiatom separations deserve further study. 
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The means used for and the feasibility of suspending antimatter in 
a vacuum depend on the physical state and properties of the antimatter. 

The means and feasibility depend on whether the antimatter is in the form 
of a plasma, gas, liquid, or solid, and whether it is magnetized, magnetic, 
paramagnetic, diamagnetic, electrically conductive, or electrically 
polarizable. Means exist for the complete or partial suspension of anti- 
matter in any of the above states and with any of the above properties, but 
a study should be made to ascertain the practical feasibility of suspension 
under the various possible combinations of states and properties. 

In addition it would be worthwhile to study the effects of annihilation 
rates on various means to transfer the antimatter from storage to the 
combustion chamber, and on the design of the combustion chamber. 
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